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I. INTRODUCTION 


In this report we consider indirect probes of supersymmetry (SUSY). Following our charge, we will review 
the current experimental status and discuss possible levels of improvements in various measurements sensitive 
to supersymmetry through either virtual or astrophysics effects. We mention the upcoming experiments which 
are likely to achieve such precision, and outline which theoretical models and ideas can be tested by those 
experiments. 

Since it is impossible to give a detailed review of every single topic here, we have limited our discussion to a 
few representative topics which were studied either at Snowmass or since then. A large fraction of this report 
is based on the individual written contributions 011S § I to our subgroup as well as talks presented at 
Snowmass. 

In Section || we discuss the implications of the recent — 2 measurement (and its possible improvements) 
for supersymmetry. In Section n we review searches for CP violation and Section deals with lepton-flavor 
violation (LEV). Section^ is devoted to the future of B-physics. In Section VI (VII) we discuss di rect ( indirect) 
searches for supersymmetric dark matter. We summarize and present our conclusions in Section VIII. 


II. ANOMALOUS MAGNETIC MOMENT OF THE MUON 

Undoubtedly, among the most exciting news of the year was the announcement of the new measurement of 
the muon anomaly at Brookhaven 0. The muon anomalous magnetic moment was reported^ to differ from 
the Standard Model (SM) prediction by 2.Qa 

a"'’ - a™ = (43 ± 16) x 10-^° , (1) 

which is about three times larger than the Standard Model’s electroweak contribution . Deviations of roughly 
this order are expected in many models motivated by attempts to understand electroweak symmetry breaking. 
A supersymmetric interpretation is particularly attractive, since supersymmetry naturally provides electroweak 
scale contributions that are easily enhanced (by large tan /3) to produce deviations of the required magnitude. 
In addition, is both flavor- and CP -conserving. Thus, while the impact of supersymmetry on other low 
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FIG. 1: Allowed values of Mlosp, the mass of the lightest observable supersymmetric particle (LOSP), and 
from a scan of parameter space with Mi = M 2 / 2 , = 0, and tan/3 = 50. Green crosses (red circles) have smuons 

(charginos/neutralinos) as the LOSP. A stable lightest supersymmetric particle (LSP) is assumed. Relaxing the relation 
Ml = M 2/2 leads to the solid envelope curve, and further allowing arbitrary left-right smuon mixing (large A^) leads 
to the dashed curve. The envelope contours scale linearly with tan/3. The \a (dark shaded, blue) and 2a (light shaded, 
yellow) allowed ranges are shown, and the discovery reaches of linear colliders with = 500 GeV and 1 TeV are 

given by the vertical blue lines. (From Ref. |^.) 


energy observables (see sections HI and IV) can be highly suppressed by scalar degeneracy or small CP-violating 
phases, supersymmetric contributions to are generic. 

Supersymmetric contributions to have been explored for many years 0 , ig, 1 ^, 1 ^, n n n, n, 
^ . Following the recent a„, result, the implications for supersymmetry have been considered in numerous 

1 ^, 0, @ n, 0 n, 13 n, n, 13 1 ^, n. The 
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studies 1 ^, ^ 1 ^, ^ 

most significant consequence is that at least two superpartners cannot decouple, if supersymmetry is to explain 
the deviation, and one of these must be charged and so observable at colliders. Non-vanishing thus imply 


upper bounds on the mass Mlosp of the lightest observable superpartner. 

Fig. 0 shows the results from a series of high statistics scans in the relevant supersymmetric parameter space, 
consistent with slepton-flavor conservation: Mi, M 2 , /r, tan/3, and A^. (For more details, see 

Ref. [^.) The points are obtained by assuming gaugino mass unification Mi = M 2 / 2 , fixing A^ = 0 and 
tan/3 = 50, and scanning over the remaining parameters up to 2.5 TeV. Collider bounds from supersymmetry 
searches are enforced and a neutral LSP is assumed. Relaxing the gaugino unification assumption leads to 
possibilities bounded by the solid curve. Finally, allowing any A^ in the interval [—100 TeV, 100 TeV] extends 
the envelope curve to the dashed contour of Fig. ^ The envelope contours scale linearly with tan /3 to excellent 
approximation. 

From Fig. we see that the measured deviation in is in the range accessible to supersymmetric theories 
and is easily explained by supersymmetric effects. The case where the LSP decays visibly in collider detectors 
yields even lower bounds, and is examined in Ref. p^ . 

Such model-independent upper bounds have many implications. They improve the prospects for observation 
of weakly-interacting superpartners at the Tevatron and LHC. They also impact linear colliders: for reference, 
the discovery reach of linear colliders with ^/s = 500 GeV and 1 TeV are given in Fig. |^. In this highly model- 
independent framework, an observable supersymmetry signal is very probable at a 1.2 TeV linear collider; the 
expected improvements in measurements may significantly strengthen such conclusions. (The ultimate target 
of the BNL experiment is an experimental error of 4 x 10“^°.) Finally, these bounds provide fresh impetus for 
searches for lepton-flavor violation, which is also mediated by sleptons and charginos/neutralinos |^ . 

Turning now to specific models, we first consider the framework of minimal supergravity, which is completely 
specified by four continuous parameters and one binary choice: toq, M 1 / 2 , Aq, tan/3, and sign(/i). The 
first three are the universal scalar, gaugino, and trilinear coupling masses at the grand unified theory scale 
Mgut ^ 2 X IQi® GeV. 

In minimal supergravity (mSUGRA), sign(a®^®'^) = sign(/rMi, 2 ), so the result prefers a particular sign 
of /r relative to the gaugino masses. As is well-known, however, the sign of fj, also enters in the supersymmetric 
contributions to B ^ Xg^. Current constraints on 33 ^ require /rMa > 0 if tan /3 is large (here M 3 is the 
gluino mass parameter). Gaugino mass unification implies Mi^ 2 M 3 > 0, therefore, a large discrepancy in is 
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only possible for > 0, in accord with the new measurement. Minimal su perg ravity, and gaugino unified 

models, in general, are generally consistent with the BNL measurement 0 §, II n, §, n, 11, §, a. 

In contrast, the minimal model of anomaly-mediated supersymmetry breaking seems to be disfavored. One 
of the most striking predictions of anomaly mediation is that the gaugino masses are proportional to the 
corresponding beta function coefficients, and so Mi^ 2 -M 3 < 0. Anomaly-mediation, therefore, most naturally 
predicts < 0 [Hj HIj in contrast to the observed deviation. The dependence of this argument on the 

characteristic gaugino mass relations of anomaly mediation suggests that similar conclusions will remain valid 
beyond the minimal model. 

In summary, the recently reported deviation in is easily accommodated in supersymmetric models. Its 
value provides model-independent upper bounds on masses of observable superpartners and already discriminates 
between well-motivated models. 


III. CP VIOLATION 


CP violation is among the least understood phenomena in the Standard Model. At present, CP violation is 
observed in only a small number of processes, such as in Kaon and B-meson mixing and decays. This can be 
accommodated through the single phase that appears in the CKM matrix. However, the CKM phase is not the 
only parameter that can lead to CP violation in the SM. The QCD 9 term, 
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327r2 


CG 


( 2 ) 


where G is the gluon field strength and G is its dual, also leads to CP violation, and indeed there are very strong 
constraints on this from limits on the electric dipole moment of the neutron (and mercury). Furthermore, CP 
violation is an essential ingredient of almost all attempts to explain the matter-antimatter asymmetry of the 
universe , yet the amount of CP violation present in the CKM matrix is insufficient to explain the observed 
, 1^ H . Hence, searches for CP violation beyond the CKM matrix are an important probe 
into physics beyond the Standard Model. 

Electric dipole moments (EDMs) violate both parity (P) and time reversal (T) invariance. If CPT is assumed 
to be an unbroken symmetry, a permanent EDM is, then, a signature of CP violation. A non-vanishing 
permanent EDM has not been measured for any of the known elementary particles. In the Standard Model, 
EDMs are generated only at the multi-loop level and are predicted to be many orders of magnitude below 
the sensitivity of foreseeable experiments ||^, H|. A non-vanishing EDM, therefore, would be unambiguous 
evidence for CP violation beyond the CKM matrix, and searches for permanent EDMs of fundamental particles 
are powerful probes of extensions of the Standard Model. In fact, current EDM bounds are already some of 
the most stringent constraints on new physics, and they are highly complementary to many other low energy 
constraints, since they require CP violation, but not flavor violation. In this Section we review the experimental 
prospects for the EDMs of various systems, and the implications for supersymmetry. 
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A. Experimental limits on EDMs 

Eirst, we summarize the current experimental bounds on EDMs and briefly mention future prospects. For 
the electron EDM, the current bound is |H 

de < 4 X 10“^^ e cm , (3) 

although a new result was recently informally announced, de < 1.5 x 10“^^ e cm [|^. It is unlikely significant 
further improvements can be expected due to stray magnetic fields. A new approach for the electron is a YbE 
molecule method [|6| that could allow getting down to the 10“^° level within a decade |H. The best current 
bound with this method is about an order of magnitude weaker than the bound quoted above. 

Eor the neutron EDM, the current bound is 

d„ < 6.3 X 10“^® e cm , (4) 

and the expectation is to strengthen the limit down to 10“^® by 2004 (H. By 2012 the goal is to reach the 
10“^®e cm level using a cryogenic apparatus. 

Eor i®^Hg atom EDM (Schiff moment), the current bound is [|^ 

dng < 2.1 X 10“^® e cm , (5) 

from a 2001 analysis. There may be hope to improve this by another factor of two, although the systematic 
errors will need to be carefully investigated. 
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B. Supersymmetric effects 


In the most general flavor non-preserving MSSM, there are over 40 new complex phases [Q. New complex 
phases arise, for example, in the Higgs mixing mass /r, as well as in the soft SUSY-breaking terms in the 
Lagrangian: the trilinear scalar mixing masses, the bilinear Higgs mixing parameter, and the gaugino mass 
parameters. Not all of these phases are physical, however, and many may be removed by field redefinitions. 

In a wide range of supersymmetric models, these CP phases are expected to be 0{1) |^. Phases of this 
size lead to an EDM of the electron, neutron, and mercury atom that are significantly larger that experimental 
bounds. There are several possible ways to avoid these experimental constraints, with varying degrees of 
sim plic ity and naturalness. One solution is to simply take the phases to be small, of order 10“^ is sufficient 
Another possibility is to allow arbitrary phases, but push (at least some of) the sparticle masses to the 
multi-TeV region which suppresses the supersymmetric contribution Embedding supersymmetry 

in a left-right symmetric framework can also suppress the phases j^. Finally, it is possible that large phases are 
permitted due to cancellations that conspire to render the SUSY contributions to EDMs below the experimental 
limits 79, 

The dominant contributions to the lepton EDMs arise from the one-loop chargino and one-loop neutralino 
graphs. For the neutron EDM, important contributions also arise from one-loop gluino graphs and two-loop 
stop-top and sbottom-bottom graphs. The operators that contribute are the electric dipole operator 


the chromoelectric dipole operator 


and the purely gluonic dim -6 operator 






( 6 ) 

( 7 ) 




( 8 ) 


In extracting the effects of the chromoelectric and the purely gluonic operators, one can use naive dimensional 
analysis to relate | 8 ^ 


dg = j (C = -— 




( 9 ) 


where 77 *^ ~ 77 *^ ~ 3.4 and M = 1.19 GeV is the chiral symmetry breaking scale. The neutron EDM d„ is 
estimated using the S't/(6) quark model Q d„ = (|dd — jd„). (For an update using QCD sum rules, see [S^.) 

The EDM of ^®®Hg arises from the T-odd nucleon-nucleon interaction in supersymmetry [Q, induced mainly 
from the color operators with light quarks. This interaction gives rise to an EDM of the mercury atom by 
inducing the Schiff moment of the mercury nucleus. The QCD uncertainties related to this calculation are 
actually smaller than for the case of the neutron EDM. This interaction can be calculated in terms of the 
MSSM phases, with the result that there are somewhat better limits using mercury EDM than the electron 
EDM 

In fact, using the combination of the electron, neutron, and mercury EDM constraints severely limits the size 
of phases in the MSSM. Although this cannot be done in general (due to the large number of parameters), some 
information can be gleaned from a more restricted, 15-parameter MSSM. If we do not require gaugino mass 
unification at the scale where the gauge couplings intersect, then there are two independent phases in the gaug- 
ino masses since one can be made real by a U{l)ii rotation. The 15-parameter MSSM considered by Ref. ||84|] 
consists of the phases in Mi, M3, A^, A„, Ag, At, and fi, as well as the real quantities tan/3, the gaugino masses, 
a common scalar trilinear coupling, |/r|, and the sfermions masses These parameters were sampled 

using a Monte Carlo scan of nearly 10® sets of parameters, and the remaining solutions satisfying all collider 
bounds plus the three EDM constraints were found. The results are illustrated in Fig. ||, in which the bounds 
on the phases of /i, A/i, and M 3 are shown. Notice that the constraint on arg(/j,) is strengthened significantly 
as tan/3 is increased. Furthermore, while there are lightly populated diagonal bands in [arg(/r),arg(Mi_ 2 )] space 
representing cancellations between diagrams, these regions suffer from larger fine-tuning and/or are (border¬ 
line) excluded by the Higgs mass limit. Clearly the parameter space is rather strongly constrained when all 
experimental constraints are applied and some lower bound on fine-tuning is imposed. 

Finally, it is interesting to note that the operator giving rise to the electron EDM is similar to the operator 
for g — 2, and this similarity can be utilized to find relations between these phenomena. In particular, assuming 
the BNL discrepancy is explained by supersymmetry, the phase of the electric dipole operator of the electron 
can be shown to be less than 2 x 10 “^ 
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FIG. 2: Parameter sets in the 15-parameter MSSM satisfying the experimental limits on the electron, neutron, and 
mercury EDMs. The constraints on 6^ = arg(/r) are the strongest, in comparison with = arg(Mi, 3 ). Open circles 
suffer from larger fine-tuning (< 0.01), defined by the maximum allowed variation of the input parameters such that the 
point survives the cuts. Lightly shaded (open or filled) circles have rrih <113 GeV. (From Ref. [^^.) 


C. EDM of the Muon 

The field of precision muon physics will be transformed in the next few years. The EDM of the muon is, 
therefore, of special interest. A new BNL experiment has been proposed to measure the muon EDM at the 
level of 


dfj, ~ 10 e cm , (10) 

more than five orders of magnitude below the current bound 

= (3.7 ± 3.4) X 10-1® e cm , (11) 

and even higher precision might be attainable at a future neutrino factory complex 

The interest in the muon EDM is further heightened by the recent measurement ^ of the muon magnetic- 
dipole moment (MDM) |^. The EDM and MDM arise from similar operators, and this tentative evidence 
for a non-Standard Model contribution to also motivates the search for the muon EDM 1^. In fact, the 
deviation of Eq. ( 0 ) may be partially, or even entirely attributed to a muon EDM ! This is because in modern 
experiments the muon MDM is deduced by measuring (the magnitude of) the muon spin precession frequency 
in a perpendicular and uniform magnetic field. However, the spin precession frequency receives contributions 
from both the MDM and the EDM. For a muon traveling with velocity /3 perpendicular to both a magnetic 
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FIG. 3: Regions in the (a”^, plane that are consistent with the observed | u; a| at the la and 2a levels. The current 
la and 2a bounds on are also shown. (From Ref. Q.) 


field S and an electric field E , the anomalous spin precession vector is 

a = -a^— B -dj^ 13 X B - E - ^ ^ ) f3 x E . (12) 

trifj, n n m^c \Y - 1 J 

In recent experiments, the last term is removed by running at the ‘magic’ 7 « 29.3, and the third term is 
negligible. For highly relativistic muons with | /3 | « 1, then, the anomalous precession frequency is found from 



where NP (SM) denotes new physics (Standard Model) contributions, and assuming d^^ ^ 

We see that the effect (|^) can also be due to a combination of new physics MDM and EDM contributions. 
Fig. H shows the regions in the (aj^^, d^^) plane that are consistent with the observed deviation in | a; a|- The 
current Icr and 2cr upper bounds on d^^ are also given. It is evident from the figure that a large fraction 
of the region allowed by both the current measurement (|^ and the df^ bound ([ll] ) is already within the 
sensitivity of phase I of the newly proposed experiment (with sensitivity ^ 10“^^ e cm). 

The proposed dedicated muon EDM experiment will use a different setup, by applying a constant radial 
electric field. In that case a similar EDM <-> MDM ambiguity is present Q, and can be resolved by up-down 
asymmetry measurements. 

It is useful to write the new physics contributions to the EDM and MDM operators as 

^NP _ e 2mA , = Tie A , (14) 

with A = 1^1 This defines an experimentally measurable quantity (j)cp which quantifies the amount of 

CP violation in the new physics, independently of its energy scale. Upon eliminating |^|, one finds 

NP 

df = 4.0 X 10-22 e cm ^3 tan,^cp • (15) 

The measured discrepancy in | tu ^1 then constrains 0 cp and d^^. The preferred regions of the {(j^cPjd^^) 
plane are shown in Fig. For ‘natural’ values of 4>cp ^ Ij dJjj^ is of order lO-^^ e cm. With the proposed dJjj^ 
sensitivity of (0), all of the 2(7 allowed region with ^cp > 10 ^ r yields an observable signal. 

At the same time, while this model-independent analysis indicates that natural values of ^cp prefer d^^ well 
within reach of the proposed muon EDM experiment, very large values of d^^ also require highly fine-tuned 
(pCP- For example, we see from Fig. I that values of dJjj^ ^10 20 e cm are possible only if |7r/2 — (?!)cp| ~ 10 
Furthermore, in specific supersymmetric models it is difficult to achieve values of (f^ large enough to affect the 
conventional interpretation of ([I]). For example, in supersymmetry, assuming flavor conservation and taking 
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FIG. 4: Regions of the (t/icp, plane allowed by the measured central value of | a; a| (solid) and its Icr and 2g 
preferred values (shaded). The horizontal dot-dashed line marks the proposed experimental sensitivity to The red 
horizontal solid lines denote the current Icr and 2cr bounds on (From Ref. |^.) 


extreme values of sparticle masses (~ 100 GeV) and tan [3 (tan /3 ~ 50) to maximize the effect, the largest 
possible value of is ~ 10“^ 1^. Very roughly, one therefore expects a maximal only of order 

(e?i/2m^c)a“®'’^ ~ 10“^° e cm in supersymmetry. Similar conclusions hold for specific models as well js^ . 
Simplest models relate the EDMs of the electron and muon by ‘naive scaling’: 


d 



(16) 


Given the current published bound on the electron EDM, ‘naive scaling’ limits the muon EDM as 

df_i < 9.1 X 10“^® e cm , 


(17) 


at the 90% CL, barely below the sensitivity of ( 0 )- Naive scaling must be violated if a non-vanishing is to 
be observable at the proposed experiment, and this may happen in one of three ways: 

• Departure from scalar degeneracy, i. e. generation-dependent slepton masses . 

• Departure from proportionality, i. e. the A terms do not scale with the corresponding fermion mass . 

• Flavor violation, i.e. non-vanishing flavor off-diagonal elements for the sfermion masses and the A- 
terms 

Such multitude of possibilities provides sufficient motivation and relatively good prospects for a dedicated 
muon EDM experiment. 


IV. LEPTON FLAVOR VIOLATION 

One of the most powerful probes of low energy supersymmetry are the precise measurements and limits on 
flavor-violating processes. In the squark sector, for example, the smallness of Kq Kq mixing either requires 
tiny off-diagonal squark (mass)^ elements, or pushes supersymmetry to embarrassingly high scales ^ 0(100 TeV) 
PA Similarly strong constraints are also present for certain elements in the slepton mass matrix. Here, 
we discuss a few of the strongest constraints from limits on rare p ^ e processes. 


A. Experimental status 

There are also several experimental probes that tightly constrain lepton-flavor violation. The strongest 
constraints arise from the rare /i —> e processes: p —*■ ey, p —> 3e, and p e conversion. The current bounds 











FIG. 5: One example set of chargino-sneutrino contributions to muon g — 2 and /r —> 67 in the interaction eigenstate basis 
with incoming right-handed muons. The photon (not shown) is emitted from the chargino. The chirality flip is shown 
by the x on the fermion line, while the the lepton flavor violating mass insertion is shown by the (8). (From Ref. [^.) 


on these processes are 


BR{fi 


BR{fi 67 ) < 
BR{n 3e) < 
e conversion) ^ 


1.2 X 10"“ @ 
1.0 X 10-12 @ 
2 X 10-12 @ , 


with the precise bound on fj. e conversion dependent on the particular host scattering nucleus. 

Improvements to these bounds are on the horizon. At PSI, there is a proposal to improve the limit on /r —> 67 
down to 10-1^ 1^, which is likely to get at least to the level of a few xlO-i^, since PSI has a x20 advantage 
over LAMPF in duty cycle and the planned detector is probably at least as capable as that of MEGA (the 
previous experiment). This experiment could possibly be pursued further at BNL, since a very intense muon 
beam will be built for another experiment. 

There is also a detailed proposal by MECO to improve the limit on /r ^ e conversion down to order 
IQ-i® at the BNL AGS. No other facility can compete at the moment, but toward the end of the decade the 
Japanese Joint Project might get involved. There is apparently a plan to build a dedicated cooled muon beam 
there at an early stage, which they claim would make possible an experiment at the IQ-^® level. Note that 
one can talk about such amazing sensitivity for this process because it has an excellent signature that does not 
require a coincidence and so is robust against high rates. 


B. Supersymmetric contributions 


Supersymmetric contributions to these processes dominantly arise when flavor violating contributions to the 
slepton mass matrix are present. The sizes of these flavor-violating elements of the slepton mass matrix are 
arbitrary in the MSSM. Particular models of supersymmetry breaking and mediation can in some cases predict 
the size of slepton-flavor violation. Some examples include supersymmetric unified theories (e.g., see |5^), and 


102 |. 


also exponentially suppressed contributions from sequestering [ 100 , 101 , 

The contributions to /x —> ey are particularly interesting, due to the strong resemblance between the operators 
giving rise to muon g—2 and those giving g, ej j^, This correspondence is best revealed diagrammatically. 
In particular, there is a precise correspondence between the diagrams that contribute to /x ^ ey (and r ^ /xy) 
with those that contribute to the muon anomalous magnetic moment [^ . This is illustrated for one class of 
diagrams involving charginos and sneutrinos in Fig. |^. Other than emitting an electron instead of a muon, the 
contribution has the same form except for the essential addition of a sneutrino flavor-mixing mass insertion and 
a propagator for the electron sneutrino. Using just this process, is it easy to see that the amplitudes for g — 2 
and /X ^ ey are related by 


mz 


(18) 
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in the mass-insertion approximation, assuming for this example . Here, is the off-diagonal 

element in the sneutrino (mass)^ matrix. One can systematically go through all classes of one-loop diagrams 
and all possible sparticle-mass hierarchies to find 


BR{ijL —!■ ey) ~ 10 


-4 


4.3 X 10-9 


m 


efi 




(19) 


where m = Max[m^, ] is the largest mass in the loop. So, if the deviation in the muon g — 2 measured at the 
BNL experiment is interpreted as supersymmetry, we obtain strong model-independent bounds on the sfermion 
mass mixing. 


TOg^/m^ < 2 X 10 ^ 

< 0.1 . 


( 20 ) 

( 21 ) 


(The bounds on the “left-left” and “right-right” slepton masses are essentially the same, using the approach of 
[ES.) The bound on the second- third generation mass mixing is much weaker due to the weaker limit on the 
decay BR{t ^ 7 ) < 1.1 x lO”® |1^|. These results are largely insensitive to any supersymmetric parameters. 
Indeed, for some mass ranges the bounds can be as much as an order of magnitude better than quoted above. 
The only assumption is that there are no accidental cancellations resulting from summing over diagrams, and 

even then the expectation is that the bounds are only mildly relaxed. _ 

Finally, there are exciting prospects for directly observing slepton-flavor physics at a linear collider |104| . 
Even with a limited number of slepton states, the observation (or non-observation) of slepton-flavor violation 
is expected to provide important clues to the underlying supersymmetry breaking structure. The reason that 
observable sflavor violation at a LC is possible is easy to understand. In many proposals to solve the SUSY- 
flavor problem, a high degree of degeneracy among sleptons (and squarks) is predicted. As a result, there is 
the potential for substantial mixing of flavor eigenstates. This can lead to substantial and observable sflavor 
violation. To be readily observable, it is necessary that the mass splittings between the states not be too much 
smaller than the decay widths, and that the mixing angles not be terribly small. In the case when the off- 
diagonal mass terms ar e comparable or larger than the widths, dramatic collider signatures are possible (see, 
e.g.. Ref. @ ^). 


V. B PHYSICS 


Rare decays and CP violating asymmetries provide another interesting hunting ground for SUSY-mediated 
processes. In the B system many rare decays involve b s{d) transitions, which are flavor-changing neutral- 
current (FCNC) processes that are forbidden in SM at tree level but occur at loop level. The one-loop processes 
involve gluonic, electromagnetic or weak penguin diagrams as well as box diagrams. Though suppressed in SM 
they are relatively large because of the CKM structure and the top-quark dominating the loop. However, SUSY 
processes may become competitive and interfere with those in SM. Depending on the sign of the interference 
term enhanced or depleted branching fractions are obtained. Due to the presence of new weak phases SUSY 
processes may affect CP asymmetries as well. While CP asymmetries of rare decays to flavor eigenstates are 
typically small in SM (< 1%) enhancements up to 20% are possible in SUSY models. For CP asymmetries 
of B decays to CP eigenstates, which are quite sizable in SM, SUSY processes either may enhance or deplete 
the effect. Other interesting penguin modes are B^ and B^ decays into two charged leptons, which are highly 
suppressed in SM and, therefore, bear a high sensitivity for New Physics. To establish a coherent picture of 
B decays and uncover SUSY contributions it is important to perform several high-precision measurements of 
rare-decay branching fractions and CP asymmetries. 


CLEO was the first experiment to observe and study FCNC B decays |1^|. However, the CLEO data sample 
is limited to 9.1 fb”^. At the asymmetric B factories, which started operation in 1999, BABAR and BELLE 
have recorded already data samples of 60 fb”^ and 44 fb”^, respectively. By summer 2002, BABAR expects 
100 fb”^. This will be increased to 500 fb”^ by summer 2005. If luminosity upgrades are successful as 
planned each experiment should reach 1 ab”^ by 2010. Experiments at the Tevatron (CDF, DO) will augment 
B samples. However, starting 2006 high-precision measurements are expected from BTEV at the Tevatron and 
LHCb, ATLAS and CMS at the LHC. Furthermore, there are ongoing discussions about a super B factory 

which would deliver 10 ab”^ 


per year [H 109|. Thus, high- 


operating with peak luminosities of 10 ^® cm-^s”^, 
statistics B samples should be available in the future. In the following we summarize expectations regarding 
radiative penguin decays |^, CP violation in B decays |10S| and B decays into two charged leptons. Whenever 
possible, the extrapolations are based on present measurements. 
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A. Radiative Penguin Decays 


Radiative penguin decays involve electroweak penguin loops or box diagrams. The largest decay is B ^ 
which is dominated by the magnetic-dipole operator O 7 . The SM decay rate contains the squares of the CKM 
matrix elements | Vts \ and the Wilson coefficient Ct- The latter accounts for all perturbative QCD contributions. 
Due to operator mixing an effective Wilson coefficient results. The non-perturbative contributions are absorbed 
into the hadronic matrix element of the magnetic dipole operator. Because of large model uncertainties one 
avoids the calculation of the hadronic matrix element by using the approximation that the ratio of decay rates 
of & —> 57 and b cev at the parton level is equal to that at the meson level (quark-hadron duality). SUSY 
processes yield additional contributions C"®™ and Cg®*", where the latter arises from SUSY operators that are 
equivalent to the chromomagnetic dipole operator Os- 

Th e bra nching fraction in next-to-leading order (NLO) in SM is predicted to be B{B Xgj) = (3.28±0.33) x 
10“^ ly^l- Recently, however, Gambino and Misiak argued for a differ ent c hoice of the charm-quark mass, which 
increases the branching fraction to B{B Xg'y) = (3.73±0.3) x 10“"* The present theoretical uncertainty 

of ^ 10 % is dominated by the mass ratio of the c-quark and &-quark and the choice of the scale parameter fif,. 
In an updated analysis using the full sample of 9.1 fb“^ CLEO has measured B{B Xg"f) = (3.21±0.43(sta()± 
X 10“"* 1 112 1, where the errors represent statistical, systematic, and theoretical uncertainties. 


0.27 


+0.18 


(sys)_o ,10 {th)> 


respectively. Because of the large errors this is consistent with the SM NLO prediction. Note that the signal 
region is dominated by continuum (75%) and BB (12%) backgrounds which have to be subtracted. The present 
relative statistical error is 13.4%. Assuming that measurements improve with luminosity as l/'s/U the relative 
statistical error a^s/B will be reduced to 4% (1.3%) for luminosities of 100 fb“^ (1 ab“^). In a super B factory 
one would expect asIB — 0.4% for 10 ab“^. The present relative systematic error is 8.4%. It is expected that 
with increased statistics the systematic error can be reduced substantially by using appropriate data selections 
and by improving measurements of the tracking efficiency, photon energy, photon efficiency and B counting. 
For 10 ab“^ the hope is to reach a systematic error of 1 — 2%. 

The CLEO B{B Xg^) measurement already provides a significant constraint on the SUSY parameter 
space. For example, the new physics contrib utio ns to B ^ Xgj, (7"®'“ and Cg®™, have been calculated using 
the minimal supergravity model (SUGRA) [113|. Many solutions have been generated by varying the input 
parameters within the ranges 0 < toq < 500 GeV, 50 < mi /2 < 250 GeV, —3 < Aojmo < 3 and 2 < tan/3 < 50, 
where a common scalar mass toq for squarks and sleptons, a common gaugino mass mi /2 and a common 
trilinear scalar coupling Aq is assumed in SUGRA. As usual the ratio of vacuum expectation values of the 
neutral components of the two Higgs doublets is parameterized by tan /3. The top-quark mass was kept fixed 
at mt = 175 GeV. Only solutions were retained that were not in violation with SLC/LEP constraints and 
Tevatron direct sparticle production limits. For these the ratios Rr = Cj ® ™ ( M w) ( Mw) and i?g = 
Cr^{Mw)ICl^{Mw) were determined. The results are depicted in Figure ^ |114[ . The solid bands show the 
regions allowed by the CLEO measurement. It is interesting to note that many solutions are already in conflict 
with the data. However, due to the theoretical uncertainties it will be difficult to uncover SUSY contributions 
at high luminosities, if the central value remains closely to the present result. 

The exclusive decay rate for B —> K*^ involves the hadronic matrix element of the magnetic dipole operator, 
which in general is expressed in terms of three g^-dependent form factors Ti(g^). Eor on-shell photons Tg vanishes 
and T 2 is related to Tg. Eor the determination of the form factors various techniques are used, introducing 
additional theoretical uncertainties. Recently, two NLO calculations were carried out, predicting SM branching 
fractions of B{B K*j) = (7.1^2^) x 10“® [ |115| and B{B K*j) = (7.913;q) x 10“® |llf ]■ The most precise 
branching-fraction measurements of the exclusive decays i3° ^ and B^ K*~^j have been achieved by 

BABAR. Utilizing kinematic constraints in the B rest frame provides a substantial reduction of gg-continuum 
background here. In a sample of 20.7 fb“^ BABAR measured S(i3° ^ K*^j) = (4.39 ± 0.41 ± 0.27) x 10“® in 
the final state [I17|. Due to the large theoretical errors of 35 — 40% the BABAR measurement is still 

consistent with the NLO SM predictions. Note that the combined statistical and systematic error is already more 
than a factor of three smaller than the theoretical uncertainty. The precision expected for increased luminosities 
will be comparable to that in the inclusive mode. Thus, it will be difficult to use the exclusive modes for SUSY 
discoveries, unless the theoretical errors are considerably reduced or SUSY effects are gigantic. In hadron 
collid ers B is also measurable. GDF expects to achieve a 7.6% statistical error per 2 fb“^, while BTEV 

[118| and LHCb [119| estimate a statistical error of crg/S ~ 0.6% per year of LHC running (~ 2 fb“^). 


Other interesting radiative penguin decays are the B Xg£~^£ modes, where £^ is either an 


or a 

In SM, these decays are suppressed by about two orders of magnitude with respect to corresponding 
B —> Xgj modes. The suppression by a is compensated partially by additional contributions from the Z^- 
penguin diagram and a box diagram that involves the semileptonic operators, Ogv and Oioa- Each of them 
can receive additional SUSY contributions. The branching fractions of the inclusive modes in S M in NLO ar e 
predicted to be B{B ^ XgC+e-) = (6.3tJ;g) x 10"® and B{B A,/x+/r-) = (5.7±0.8) x 10"® pOt 
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FIG. 6: Scatter plot of Rs, versus R7 for solutions obtained in the SUGRA model. The region allowed by the CLEG 
measurement lies inside the two sets of solid diagonal bands. 


The theoretical uncertainties are about 16%. These modes have not been observed so far. The lowest branching- 
fraction upper limits @90% from CLEO are about an order of magnitude above the SM predictions 1123 . Using 
SM predictions and efficiencies determined by CLEO an observation of the Xse^e~ and modes is 

expected in a sample of 100 fb“^ with statistical errors around ^ 17% and ^ 19%, respectively. This will be 
improved to ^ 5.3% 1.7%) and ^ 6% 1.9%) in samples of 1 ab“^ (10 ab“^), respectively. Unless the 

SUSY contributions lead to significant enhancements the theory errors need to be reduced at the same time 
precise measurements are obtained in order to increase the sensitivity for observing New Physics. 

Branching fractions of the exclusive modes are further suppressed. Using SM predictions from two recent 
models and their uncertainties yield the following ranges of branching fractions: B{B K£~^£~) = (4.7 — 
7.5) X 10"^ B{B K*e+e-) = (1.4 - 3.0) x 10"®, and B{B ATV+m”) = (0.9 - 2.4) x 10“® 

SUSY processes could enhance these branching fractions considerably. As an example Figure ^ depicts the 

'' B ^ K*ijX iJL~ calculated in SM, SUGRA models and minimal-insertion- 
Shown are both the pure penguin contribution and the sum of the penguin 


approach SUSY models (MIA) 125 


process and the long-distance effects, displaying constructive interference below the charmonium resonances and 
destructive interference above. The different models are characterized in terms of ratios of Wilson coefficients 
Ri = \ for i = 7,9, 10. The SM prediction is the lowest but bears large uncertainties. 

Except for an unconfirmed signal in B ^ K'^jjX seen by BELLE, none of the exclusive B K{K*)£^£~ 
modes have been observed yet. The branching fraction oi B{B = (0.99to;32lai4) ^ measured 

by BELLE |126| is barely consistent with the BABAR limit of B{B —> K£~^£~) < 0.6 x 10“® @90% CL ]117| . 
The BABAR limits of B{B K*°e+e-) < 5.0 x 10"®, and B{B < 3.6 x 10"® lie less than 

a factor of two above the SM predictions. In a sample of 100 fb“^ we expect first observation of these modes. 
The statistical errors expected at high luminosities for B —> K*^£^£~ are about a factor of two higher than 
those for the corresponding inclusive mod es. Experiments at the Tevatron and LHC will be competitive in the 
Ar*°/i+/r“ and final states |118| . 

The lepton forward-backward asymmetry Afb{s) as a function of s = is an observable that is very sensitive 
to SUSY contributions. It reveals characteristic shapes in the SM both for inclusive and exclusive final states. 
To avoid complications from the charmonium resonances one restricts the range s to masses below the 
which accounts for ~ 40% of the entire spectrum. Figure 0 shows Afbis) for the B mode |125| . 

In SM, the position Sq of Afbiso) = 0 is predicted to lie at sq = 2.881q 28 GeV^. Both, the shape and sq are 
expected to differ significantly in New Physics models. The shape is very sensitive to the sign of Rr and varies 
from model to model. Thus, a precise measurement of Afb{q^) may permit an extraction of the coefficients 
Ri. However, to achieve sufficient precision a super B factory is needed. For Measuring 18 data points below 
s = 9 GeV^ with 100 events each in the B —> Xs£'^£~ {B K*°£'^£~) modes at a super B factory (10 ab“^/y) 
requires a run period of 0.3-0.4 (0.8-1.3) years. For comparison, LHCb expects to achieve the same precision in 
the B mode in about one year. 

The modes bearing the smallest theoretical uncertainties are the B Xgvv decays, since long distance 
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FIG. 7: The dilepton invariant mass-sauare d sp ectrum (left) and the normalized forward-backward asymmetry (right) 
as a function of s = in B ^ K*fj,~ [[125f| . The solid lines denote the SM prediction. The shaded region depicts 
form-factor related uncertainties. The dotted lines correspond to a SUGRA model {R 7 = —1.2,Rq = 1.03, Bio = 1) and 
the dash-dotted lines to a MIA model {R 7 = — 0.83,B9 = 0.92, Bio = 1.61). In the spectrum both the penguin 
contribution and the distribution including long-distance effects are shown. In the Afb plot the upper and lower sets of 
curves show the difference between ^ < 0 and > 0, while the dashed curves give results for another MIA 

model (B 7 = =f0.83, Bg = 0.79, Bio = —0.38). 


effects are absent and QCD corrections are small. Here, only weak penguin and box diagrams contribute. 
However, experimentally both inclusive and exclusive modes are difficult to observe due to the two neutrinos. 
To reduce backgrounds from qq continuum and other B decays at least a partial reconstruction of the other 
B is nece ssary . _ Si nce branc hing fractions are slightly lower than those in the corresponding B 


channels 124, 127, 128, 129|, several years of running at a super B factory (at 10^® cm ^s ^) 




are needed to 

study these modes. For example, in a sample of 10 ab“^ the relative statistical error expected in the inclusive 
B Xgvv branching fraction is as IB ~ 5%. For the K*vv and Kvv final states the relative statistical errors 
are crg/H ~ 10% and asjB ~ 13%, respectively [0. 


B. CP Violation 


CP violation can arise from different effects. We focus here on (i) CP violation resulting from the interference 
between decays with and without mixing, which occurs when decay rates of a and into a CP eigenstate 
or into CP-conjugate states differ, and on (ii) CP violation in decay, resulting when the magnitudes of the 
amplitudes for a decay and for the CP-conjugate decay are different. In SM, the only weak phase besides the 
mixing phase is the phase of the CKM matrix. In SUSY-mediated processes new phases arise as well as new 
contributions to P°P° mixing. Experimentally, CP violation is present when an asymmetry between a process 
and its CP conjugate is observed. At the T(4S') the two B mesons are produced in opposite CP eigenstates. 
For such systems time-integrated CP asymmetries of P° and P° mesons decaying into a CP eigenstate or into 
CP-conjugate states vanish. To observe an effect, one needs to measure the CP asymmetry as a function of the 
time difference At between the two B decays. Measurements of such time -dependent CP asymmetries provide 
a determination of the angles a, /3 and 7 in the Unitarity Triangle |130 . For measurements of 7 also other 
methods exist. The focus here is on measurements accessible in neutral and charged B decays. 

The time-dependent CP asymmetry for a B decay into a CP eigenstate, fcp, is defined by 


Mcp(At) = 


r(P° ^ /cp)(At) - r(p° ^ fcp){At) 
r(P° ^ /cp)(At) + r(po ^ /cp)(At) 


= Cf cos(AmB^ At) -I- Sf sin(AmB^ At), (22) 


where 

























13 


Cf = 


(1-1 A P) 
(1+IAP)’ 


Sf = 2 


ImX 


(1+|AP)^ 


and A = 


Vf 


p A 


(23) 


The factor rjf indicates the CP eigenvalue of fcp, yielding pf = ±1 for CP{fcp) = ±1, q/p represents the 
mixing contribution and A/A is the amplitude ratio for ^ fcp and B^ —> fcp, respectively. CP is 
violated if | A 1 or IttiX ^ 0. Experimentally, the measurement of Acp{^t) involves three steps: (i) the 


reconstruction of CP eigenstates states such as B ^ or B 


(ii) the tagging of the h flavor 


at production point using for example the charge of a lepton or a kaon observed in the other B meson, and 
(iii) the measurement of the time difference At of the two i?-decay vertices. To account for time-resolution 
effects, the measured CP asymmetry is parameterized by the periodic function in equation ^ convoluted with 
a time-resolution function. To account for errors in the tagging procedure, the measured amplitudes C/ and Sf 
contain in addition a dilution factor D = 1 — 2w, where w represents the fraction of mistagged events. 

For B J/jjjKg and related modes one expects | A |= 1 in SM. Thus, the first term in equation ^ 
vanishes and CP asymmetries measure D ■ Sj with Sj = sin 2/3. Using a sample of 29.7 fb“^ BABAR was 
the first to observe CP violation in the B^B'^ system. The CP sample contained 803 events of which 640 
events remained after tagging and vertexing. By analyzing CP asymmetries of B^{B^) decays into J/il’Kg, 
tj;{2S)Kg, XcKg, and J/ij^Kg, CP eigen states as well as CP conjugate states, BABAR measured 

sin 2/3 = 0.59 ± 0.14(s(jj() ± 0 .05(- s„.s~i 1131]. BELLE confirmed the presence of CP violation measuring sin 2/3 = 
0.99 ± 0.14(stai) ± 0.06(sys) [ll32| . The present world average of sin 2/3 = 0.79 ± 0.1 is consistent with a value 
of (3 obtained from measurements of the three sides of the Unitarity Triangle. However, presently the errors 
are rather large, mainly because of large theoretical uncertainties in the extraction of CKM parameters from 
measurements of semileptonic branching fractions, Ck (parameterizing CP violation in mixing), and 

AmB^- The measurement of a is complicated by contributions from penguin processes. Thus, both time- 
dependent terms in equation |2^ are present and ImX = sin2ae//, where a and cteff differ by a penguin phase 
6 p. Furthermore, a-sensitive modes are b ^ u transitions that are suppressed with respect to b c transitions 
by I VubIVcb P= 0(10“^), thus requiring larger B samples than for 6 ^ c processes. So far a has not been 
measured directly. In a sample of 30.4 fb“^ BABAR has studied the time dependence of 65 tt"* 


yielding 5.. = 0.03+;^J^ ± and C.. = -0.25+^o;f, ± 0.14(,,,) 

The measurements of the three angles in the Unitarity Triangle allow various tests. First, one checks for 
consistency between /3 measured directly in CP asymmetries and /3 obtained from measurements of the sides of 
the Unitarity Triangle. Second, one compares sin 2/3 measurements obtained in different quark processes, such 
as b ^ CCS, b ccd and b sss. In SM, all measurements have to give the same result. SUSY contributions, 
however, may affect each quark process differently. Thus, it is necessary to measure sin 2/3 also in 6 —> ccd 
processes such as B ^ £){*)+£)(*)- qy B ^ J/'(/i 7 r° and in 5 ^ sss processes such as B ^ equally well as in 
B —> J/tpKg. Third, with the additional measurements of a the Unitarity Triangle is two-fold overconstrained. 
Fourth, CP asymmetries measured in i3 —> J/ipK'g and B —> tt+tt”, combined with measurement of the 
B^B^ oscillation frequency Ams^ and that of | VubIVcb \ allow a model-independent analysis to look for New 
Physics contributions. These can be parameterized in terms of a scale parameter that 

accounts for new contri butions to B^B^ mixing, and a new weak phase 9d that represents new sources of CP 
violation |134, |135|, 1361. Finally, with the measurement of 7 one can test for closure of the Unitarity Triangle: 
a -f /3 -I- 7 = TT ||13 6 || . 

In order to carry out these tests and to uncover SUSY contributions in CP asymmetries high-precision 
measurements are important. The present statistical uncertainty of sin 2/3 in BABAR is 0.14 for 29.7 fb“^. The 
precision is expected to improve as Assuming that the reconstruction efficiency and tagging performance 

remain unchanged extrapolations to 100 fb“^ and 1 ab“^ yield statistical errors of crsin 2/3 = 0.076 and (Tsin 2/3 = 
0.024, respectively. Since BELLE should achieve similar precisions, the error in the world average of sin 2/3 
is reduced by another factor of \/2. Contributions of the systematic error consists of tagging uncertainties, 
vertexing resolution, B life-time precision and background treatment. With increased statistics, the individual 
systematic errors will be reduced, so that the total systematic uncertainty should remain of the order of the 
statistical error. For an annual luminosity of 10 ab“^ at a super B factory the statistical error will be reduced 
to crsin 2/3 = 0.0076. This should be compared with expectations at the LHC. Based on Monte Carlo simulations 
the statistical precision of sin 2/3 measured in CP asymmetries of i3 ^ is estimated to be crsin 2 pitpKs) = 

0.017 for AT LAS, crsin 2 / 3 (V’^s) = 0.015 for CMS a nd g sin 2 /?(UAl,s) = 0.021 for LHCb for one year of LHC 
running [1W|. In a sample of 2fb~^ at the Tevatron ||ll 8 | precisions expected for sin2/3 are crsin 2 / 3 (V'-f^s) = 0.05 
for CDF, a sin 2 pii’Ks) = 0.04 for DO, and asin 20 {'<pKs) = 0.025 for BTEV. In the next 5-10 years, when the 
high luminosities will be achieved, the measurements of the sides of the Unitarity Triangle will also improve. 
Particularly, the theoretical uncertainties associated with the side measurements are expected to improve by a 
factor of 2-4 |10£]. 
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FIG. 8: Uncertainty in extracting the angle Sp = a^ff — a in degrees versus the branching fraction for 

samples of 90 fb“^, 500 fb“^, and 10 ab“^. 


For a comparison of sin 2/3 measurements in different quark processes precisions expected for b —>■ sss and 
b ccd processes are not s uffic ient in present asymmetric B factories. In a sample of 22 fb“^ BABAR observes 
11 B ^ ipKg events |137, 13S|. Using the results of the BABAR sin2/3 measurement and the observed 4>Kg 


yield gives statistical-error estimates of (Tsin20{4‘Ks) = 0.56 for 100 fb“^, asin 2 i 3 {(l>Ks) = 0.18 for 1 ab“^, and 
crsin 2 i 3 {(t>Ks) = 0.056 for 10 ab“^. At these levels of precision syste matic errors are less important. In the 
modes B BABAR observes 13 events in a sample of 23 fb“^ |109| . Following the same extrapolation 

procedure as for (j)Kg the precision expected for sin2/3 measurements in is crsin20{J/'4’^^) = 0-52 in 

100 fb“^, crsin 2 / 3 ('^/^'^°) = 0.16 in 1 ab“^ and crsin2/3(>^/^^°) = 0.052 in 10 ab“^. 

Measurements of sin 2a also require high statistics, since branching fractions for b u processes are small 
and are affected by competing penguin amplitudes. Extrapolating the present BABAR B tt+tt" results to 
high luminosities yields crc„„ = 0.26, as„„ = 0.32 for 100 fb“^, = 0.09, = 0.1 for 1 ab“^, and 

crc^^ = 0.026, = 0.032 for 10 ab“^. To extract a from one needs to measure the penguin p hase dp. 

This is achieved by exploiting isospin relations among the amplitudes of 33 —> tttt and B tttt decays [139|. In 
the absence of electroweak-penguin amplitudes the isospin relations form two triangles (one for B and one for B 
decays) with a common amplitude base A(i3+ ^ = A{B~ 7 r“ 7 r°). The angle between the two triangles 

is 25p. The presence of electroweak-penguin amplitudes introduces a small correction. The real challenge in 
this analysis, however, is the measurement of B{B^ and S(i3° ^ since branching fractions 

are rather small and the final state is affected by a larra gg-continuum background. The determination 

of 5p as a function of jB(B^ —!■ 7 r° 7 r°) is depicted in Figure g for three different luminosities. In a 10 ab“^ 
sample the uncertainty expected for 5p is asp = 5°, which is about a factor of three larger than the uncertainty 

For 2 fb“^ BTEV estimates a statistical error of \ a: 


expected for aeff from For 2 fb ^ BTEV estimates a statistical error of = 0.024 [11B|. For 

comparison, LHCb quotes errors of = 0.07 and = 0.09 per year of LHC running. Using a Dalitz plot 
analysis in the B ^ pir channel is another promising method to extract a with high precision both at a super 
B factory and in future experiments at hadron machines. For example, LHCb quotes an expected statistical 
error of cra{pTr) =3 — 5° per year of LHC running |119|. 

One promising method of measuring 7 at the T(4S') is based on studying B decays |14C]. Here 

the different CKM structure in 6 ^ cus and b —>■ ucs processes is exploited, since the former decay has no weak 
phase while the latter involves 7 . For the extraction of 7 at least two flavor eigenstates, {K~Tr~^, A'“ 7 r+ 7 r°, 
or iF^Tr+Tr+TT") need to be reconstructed, since besides 7 the B~ branching fraction and two 

strong phases need to be measured. Note, that this procedure leads to an eightfold ambiguity in the value of 7 . 
To gain sufficient statistics actually all possible DK, D*K, DK* and D*K* modes have to be combined. In a 
sample of 600 fb“^ the statistical uncertainty in 7 is estimated to range around < 7 .^ ~ 5 — 10° depending on the 
exact value of 7 . Thus, for a range of two-three standard deviations each solution covers most of the 7 region. 
In a sample of 10 ab“^, however, this is different, as the statistical uncertainty is expected to be reduced to 
CT-), ^ 1 — 2.5°. So, again high-statistics B sam ples are required. In 2 fb“^ BTEV expects to measure 7 with an 
error of about 7°, using 33° ^ D ~decays [118|. 


Another method of measuring 7 in 33 decays involves CP asymmetries in or D^*^p modes |4l[ |^, |l4| 
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Here, the combination of CKM angles 2/3 + 7 is measured, which results from an interference between the decay 
b cud and the j?°-mi xed process b ucd, where the latter depends on the angle 7 . BABAR recently 
studied the Zl*+7r“ mode [109 , which has a large branching fraction and can be partially reconstructed with 
low backgrounds but yields a small CP asymmetry. The main disadvantage of this technique is the necessity of 
measuring the ratio r of the doubly CKM-suppressed decay b ucd to the allowed decay b cud. Extrapolating 
the BABAR study yields statistical-error estimates of crsin{ 2 f 3 +-y){D*TT) ^ 0.3 for 100 fb“^, crsin( 2 i 3 +'y)iD*'^) ^ 0.1 
for 1 ab“^, and crsi^(^ 2 / 3 +'y){D*n) ^ 0.03 for 10 ab“^. The probability that r is not measured due to statistical 
fluctuations is 30%, 10% and 3%, respectively. For comparison, LHCb expects to achieve a statistical error of 
o'(im( 2 / 3 -i- 7 ) = 0.26 per year of LHC running. 

CP asymmetries of rare decays into flavor eigenstates are also rather suited to search for SUSY- medi ated 
processes. For example in H —s- decays CP as ymm etries are expected to be small in SM (< 1%) [144| but 

they may be enhanced up to 20% in SUSY models [145|. So far all observed CP asymmetries are consistent with 
zero In the inclusive B Xsj mode CLEO measured Acp{B A^y) = (—0.079±0.108±0.022) x (1.0±0.03) 
11461, where the first error is statistical while the second and third errors represent additive and multiplicative 
systematic uncertainties, respectively. Though the CP asymmetry has been corrected for contributions from 
B — !■ Xd, a separation of these events from B —> Agy events may become necessary on an event-by-event basis 
to rule out that the CP asymmetry in H ^ A^y is canceled by that in H > A^y. In the exclusive B'^ Ar*°y 
modes BABAR observed a CP asymm etry of Acp{B Ar*°y) = —0.035 ± 0.076 ± 0.012 using the three K* 
final states, AT+tt”, and Kgir'^ jllll- Extrapolating the measured statistical error to high luminosities 

yields a^cp = 3.3% (3.5%) for 100 fb“^, a^cp = 1-0% (1-1%) for 1 ab“^ and = 0.33% (0.35%) for 
10 ab“^ in H > A^y {B A*y) modes. LHCb expects to measure Acp{B —> A*y) with a precision of 1% 
per year of running. 


C. B Decays into Two Charged Leptons 


Another interesting probe of supersymmetry is the decay Bg The SM prediction is given by 

B{Bs A^M”) = (3-7 ± 1.2) x 10“®, with the uncertainty (±25%) dominated by the decay constant fp^- The 
curre nt ex perimental bound on the branching ratio, {B), has been set during Run-I of the Tevatron, where 
CDF | |147| determined B{Bs /j,+/i“) < 2.6 x 10“® at 90% C.L. In addition to the experimental challenge, the 
almost three orders of magnitude gap between the current experimental bound and the SM prediction make s thi s 
mode an excellent laboratory for new physics. In contrast to observables which enter the unitarity triangle |148| , 
in the MSSM, the branching ratio B{Bs grows like tan®/3 |l49| , 15C , 151 , 152 , 153 |, wi th a possible 


several orders of magnitude enhancement. More interestingly, it has been very recently shown |154| that in 


the mSUGRA scenario there is a strong correlation with the muon anomalous magnetic moment {g — 2)^. An 
interpretation of the recently measured excess in {g — 2)^ in terms of mSUGRA corrections implies a substantial 
supersymmetric enhancement of the branching ratio B{Bs if {g — 2)^ exceeds the Standard Model 

prediction by {Saf_i)suSY = 4 x 10“®, B{Bs —> gAg,~) is larger by a factor of 10-100 than in the Standard Model 
and within reach of Run-H of the Tevatron. The single event sensitivity of CDF at Run-Ha is estimated to be 
1.0 X 10“®, for an integrated luminosity of 2 fb“^ ]155| . Thus if mSUGRA corrections enhance B{Bs g~) 
to e.g. 5 X 10“^, one will see 50 events in Run-Ha. Run-IIb may collect 10-20 fb“^ of integrated luminosity, 
which implies 250 -500 events in this example. 

Following Ref. [1^|, in Fig.||we show the contours of B{Bs g'^g~) (solid) and [5a^)suSY (dashed) in this 

plane, for tan/3 = 50, Aq = 0, /i > 0 and mt = 175 GeV. The shaded region is excluded through the various 
theoretical and experimental constraints. A sensitivity of B{Bs g'^g~) ^ 2 x 10“^ at CDF now corresponds 
to a sensitivity of M 1/2 ~ 280 GeV and Mq ^ 400 GeV, respectively. While GDF is not able to see squark masses 
directly up to 0.7 TeV (corresponding to M 1/2 = Mq ~ 300 GeV), it will nevertheless be able to prepare the 
ground for LHC by observing the Bg —> g^g~ mode. Even better, after 10 fb~^ CDF will probe M 1/2 ~ 450 GeV 
and Mq ^ 600 GeV (for tan /3 = 50) which in mSUGRA corresponds to masses for the heaviest superpartners of 
1 TeV. We conclude the discussion of Fig.^ wit h th e prediction of the light Higgs boson mass Mh (dot-dashed 
line) for tan/3 = 50 in the mSUGRA scenario [156]. Any measurement of B{B, n+,i- 
useful upper bound on Mh- The simultaneous information of B{Bg g 
of the (Mi/ 2 , Mo)-plane. A Higgs mass around 115.6 GeV results in 10 
would most likely be measured before the Higgs boson is discovered. 


C g ) and 5a 


B{Bg 


g^ g^) by itself implies a 
^ hxes Mh in most regions 
^ g'^g~) ~ 3 X 10“^ which 

















16 


> 

OJ 

C3 


2000 


1500 


1000 


500 


tan/?=50, Aq = 0, /x> 0, m^ = 175 GeV 



500 750 1000 

Mi /2 [GeV] 


1250 1500 


FIG. 9: Contour plots of the B(Bs —> /i^/i ) (solid) and on {5a^)susY (dashed) in the (Mo, Mi/ 2 )-plane for mSUGRA 
parameter values as shown. Contours of the light Higgs boson mass (dot-dash line) are also shown. From Ref. |154|. 


VI. DIRECT DARK MATTER SEARCHES 


Neutralinos in supersymmetry are well-motivated candidates to provide much or all of the non-baryonic dark 
matter. In many models, the neutralino is the LSP and is effectively stable. In general, the lightest neutralino 
is a mixture of the superpartners of Higgs and electroweak gauge bosons. Remarkably, detailed calculations of 
the thermal relic density of neutralinos have shown that neutralinos may indeed account for most of the missing 
mass of the universe. 

Typically, there can be several different cosmologically preferred regions in parameter space. For ex¬ 
ample, in mi nima l supergra vity, there are four: a ‘bulk’ region at relatively low mo and mi/ 2 , a ‘focus- 


point’ region [157 


158, 159 at relatively large mg, a co-annihilation ‘tail’ extending out to relatively large 


mi /2 [160, 161, 162|, and a possible ‘funnel’ bet ween the focus-point and co-annihilation regions due to rapid 
annihilation via direct-channel Higgs boson poles |163, 1641. A set of benchmark supersymmetric model param¬ 


eter choices was proposed [165] just prior to Snowmass 2001, with points representing each of those regions. In 
this and the next section, we shall use this set of benchmarks to illustrate our discussion of the various signals 
of supersymmetric dark matter |166|. 

Neutralinos are very weakly interacting and they pass through collider detectors without leaving a trace. 
Therefore, it is practically impossible to observe them in collider experiments directly. Existing bounds on 
neutralinos must rely on model-dependent correlations between their properties and those of other supersym¬ 
metric particles. However, if neutralinos make up a significant portion of the halo dark matter, many additional 
avenues for their detection open up. They may deposit energy as they scatter off nuclei in terrestrial, usually 
sub-terrestrial, detectors. 

Searches for this signal in low background detectors have been underway for around 20 years. The first 
decade was dominated by conventional HPGe (and Si) semiconductor detectors. The design of these detectors 
was to some extent “off-the-shelf” and progress was achieved, for the most part, by improving the radioactive 
backgrounds around the detectors. In the mid-90’s results from Nal scintillator detectors became competitive. 
They were able to employ pulse shape discrimination to make statistical distinctions between populations of 
electron recoil events, and nuclear recoil events. In principle, the intrinsic background of the detector and 
environment were no longer the limiting factors, since with sufficient exposure time and target mass the limits 
could be driven down. However, the relatively poor quality of the Nal discrimination meant that systematic 
effects rapidly dominate, halting any further improvement with exposure. The Nal detector technology could 
also be described as off-the-shelf, however, the low background and high light yield housing systems were very 
definitely novel. At the end of 90’s we finally saw results, from new cryogenic (<1 K) detector technology that 
had been developed specifically for direct detection, take the lead in terms of sensitivity. 

Current sensitivities are now at a WIMP-nucleon normalized cross-section of 4 x 10“® pb, which for a Ge 
target with a detection threshold of 10 keVr is an integrated event rate of 1.0 kg“^ d“^. In this discussion 
keVr will be used to indicate actual recoil energy of an event, whereas keVee (electron equivalent), will be 
used to indicate the visible energy of a nuclear recoil event based on a energy scale for electron recoil events. 
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For example, in Nal scintillator, an iodine nuclear recoil of 22 keVr, generates the same light output as a 2 
keVee electron event. To first approximation the differential recoil spectrum for WIMP nucleon events has a 
characteristic, but far from unique, exponential decay shape. This text is not the right place to d iscu s s th e 
details of what determines the spectrum, or the details of the cross-section normalization. See instead |167, 


168 


We should touch on the DAMA experiment, which has reported a 4a observation of annual modulation in the 
count rate of events in their lowest energy bins (2-6 keVee , wh i ch c orresponds to 22-66 keVr for nuclear recoil 
events scattering on iodine) over a time spanning 4 years |169|, 1701. The modulation has a phase (maximum 


in June) and a period (1 year) which would be consistent with the fluctuations in the observation arising from 
a contribution from a WIMP recoil spectrum which will be modulated by changes in the Earth’s velocity in 
the galaxy. (Baseline model assumes an isothermal, non-co-rotating WIMP population.) It is enormously 
important that this signal be studied and the source clearly determi ned . The DAMA experiment will check its 
own results with two new years of data already available for analysis |171], and further data from running for the 
seventh year just beginning. It could be argued that a substantial overhaul in the data taking strategy should 
occur before further data taking, since the significance of the modulation effect is really systematics, rather 
than statistics limited, at this stage. The importance of a “beam-off” component of the data taking cannot be 
emphasized enough to establish the long term stability of data in the lowest energy bins when looking for few 
% fluctuation. While one cannot turn the WIMP wind off to order, a revised acquisition strategy that included 
direct calibration of the stability of the acceptance of the 2-3, 3-4, 4-5 keVee bins (in which all of the WIMP 
signal is expected) would help demonstrate that the modulation signal is not an artifact. Multiple scattering 
events should be stored, rather than vetoed since they could provide a non-WIMP calibration in real time. 
Additionally, on-line light pulses, or daily gamma calibrations, could be used to demonstrate the stability of 
the count rate in the bins if even higher statistics are required. The DAMA Collaboration is now constructing 
LIBRA, which will be a larger 250 kg array which by virtue of employing lower background Nal should be able 
to make a more sensitive check of the annual modulation signal. Using a different experimental setup will also 
help address the systematics question, although its operation at the same location, by the same collaboration 
doesn’t allow all possible systematics to be checked. To this end the Boulby DM Collaboration may be able to 
carry through on the intention of running around 50 kg of Nal to cross-check the results directly with the same 
target material. Results from this Boulby program ar e un likely to be availabl e bef ore 2004, however. 

The currently reported results from the CDMS I ]172|] and EDELWEISS |173| experiments can be viewed 
as being inconsistent with the size of the annual modulation signal seen in DAMA. This is most readily seen 
by converting the published modulation amplitude of DAMA, into a cross-section on nucleon (1.4 x 10“^ pb), 
assuming scalar WIMP nucleon couplings, and then for comparison with CDMS I, to an estimated number of 
WIMP events (40) that would be observed in the 10.6 kg-live days exposure of Ge. This experiment saw 13 
single scatter nuclear recoil events in Ge during this period which is clearly inconsistent with 40. In fact, the 
observation of an additional 4 multiple scattering nuclear recoil events during the same run, can be attributed 
unambiguously to neutrons. This permits a modest reduction of the (90% CL) upper limit in CDMS I single 
scatter nuclear recoil data to 8 WIMP events. On this basis, CDMS I and DAMA are inconsistent at 99.98%. 
The Edelweiss experiment has a similar WIMP sensitivity using detectors that are like the ones used in CDMS 1. 
The data reported to date is for a smaller exposure of only 3 kg-days, however, it appears free of neutrons above 
a 25 keV analysis threshold due to the much deeper site location. Further data from longer exposure is eagerly 
awaited. 

A number of avenues exist for reducing this incompatibility between the current results of the nuclear recoil 
discriminating Ge detectors and DAMA Nal annual modulation. Firstly, the DAMA collaboration reports an al¬ 
lowed region that combines not just the positive signal from the annual modulation amplitude, but also negative 
results/upper limits from other Nal derived data. This has the effect of reducing the headline expectation for 
the WIMP cross section by approximately a factor 3. This makes the disagreement between the experiments not 
statistically signihcant, but does beg the question whether the various DAMA experimental results themselves 
are mutually compatible. Alternative routes for explaining the apparently contradictory results have looked at 
the nature of WIMPs themselves. Comparison of Na, I and Ge recoil rates requires the assumption of a model. 
As has been said before the primary analysis of the experimental data is developed assuming scalar, or spin 
independent interactions. The interaction rates of the spin independent process often dominates because of the 
enhancement by the coherent scattering across multiple nucleons. However, in the absence of spin independent 
couplings, the relative spin dependent couplings will be stronger for an Nal target, in which both elements are 
monoisotopic with odd proton spin, compared to Ge, which has only 1 odd neutron spin isotope, with an 8% 
abundance in the natural element. Existing indirect WIMP detection experiments, are able to (partially) test 
the spin dependent hypothesis, ruling out some, but not all possible solutions (see next section). Another en¬ 
tertaining model is to suggest that the WIMPs themselves inelastically scatter from nuclei (i.e. WIMPs require 
internal transition to scatter, between low energy non-degenerate states.) Tuning of the excitation energy can 
ensure that a heavier nuclei (such as I) may be able to scatter WIMPs at a significantly enhanced rate relative 
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to Ge [|174 |. 

If we now look forward at some of the predicted goals of a few experiments over the next decade, it is 
immediately apparent that the forecast rate of progress appears to be rapidly accelerating when compared to 
historical progress. The question is whether this is simply a “triumph of hope over expectation”, or represents a 
genuinely improved rate of progress that stems from applying detector technologies (2-phase Xe, semiconductor 
and scintillator cryogenic detectors, naked HPGe) that were “birthed” with this specific application in mind. 
We believe that the optimism is justified. At present, there are a number of experiments under construction 
(or that will be shortly) that will test for WIMP interactions in the range 1.0 kg“^ d“^ - 0.01 kg“^ d“^. The 
CDMS II, EDELWEISS II and CRESST II which all employ low energy threshold (< 10 keV), nuclear recoil 
discriminating, cryogenic detectors will be deploying « 10 kg of target. 

There is also significant interest in liquid Xe for WIMP detection. In principle, Xe is well placed, having a 
photon yield similar to that of Nal, while potentially being capable of creating a radioactively cleaner target 
(especially if 85Kr is removed during isotopic enrichment stage). Nuclear recoil versus gamma discrimination is 
also possible, either by pulse shape analysis, which gives a relatively weak separation, and using a simulatneous 
measurement of the photon and electron-ion charge yield in Xe which is a much stronger method of discrimina¬ 
tion. The Rome group |169 achieved early operation of Xe deep underground. A new program has been started 
by a Japanese group in Kamioka |P^|. There is a major effort on Xe by the Boulby DM Collaboration | 176 | 
which is constructing a series of liquid Xe experiments at Boulby mine in the UK. ZEPLIN I, now running at 
Boulby, is based on pulse shape discrimination. ZEPLIN II makes use of simultaneous collection of scintillation 
and charge to achieve factors of 10-100 improved sensitivity, and ZEPLIN III incorporates a high E-yield in the 
liquid to enhance the recoil signal. A new US group, XENON, has also recently been proposed, centered at 
Columbia that would also use high field operation of two phase liquid-gas Xe. 

A protot ype 1 m^ low pressure gas (CS2) detector, DRIFT, has also just started low background operation 
at Boulby |177|. 


The GENIUS [178, 179| and MAJORANA experiments have proposed to improve limits in WIMP searches 
by significant reduction in radioactive background levels, and exploitation of active Ge self-shielding. Nuclear 
recoil discrimination is not available for low energy events in conventional HP Ge detectors. It is worth noting 
that if these experiments can reduce low energy backgrounds by a factor 3000 from current levels (as projected 
by GENIUS for the complete 14 m liquid nitrogen shield) the limiting background becomes events from electron 
scattering of pp solar neutrino flux. 

In order to move beyond a sensitivity level of 2 x 10“® pb targets of 100-1000 kg (equivalent Ge) with very 
good nuclear recoil discrimination (rejecting >99.99% of electron events) will be necessary. A limiting sensitivity 
for realistic detecto r arr ays, based on existing cryogenic detectors (such as CryoArray or liquid Xe (such 

as ZEPLIN-MAX [181| will be at an event level of 1 event per 100 kg per 1 year which is equivalent to a 
WIMP-nucleon cross-section of 10“^° pb. 

Fig. |l^ shows the spin-independent cross-section for neutrali no-p roton scattering for the benchmark points 
of [165|, using two different co des: N eutdriver [168| and SSARD [182|. (Experiments sensitive to spin-dependent 
scattering have inferior reach |166[ .) One finds reasonable agreement, with th e largest differences arising for 


points D and K, where the cross-section is abnormally small due to cancellations [1183, 184[. For any given tan/3, 


the cancellations occur only for a specific limited range in the neutralino mass. Unfortunately, points D and K 
fall exactly into this category. 

Fig. 1^ also shows the projected sensitivities for a number of experiments listed in the figure caption. It 
should be emphasized that the CDMS II experimental reach, based on a modest target mass of 5 kg of Ge, will 
test a significant subset of SUSY models, and comes close to testing 4 of the benchmarks I, B, E and L. These 
could be reached by the larger (100 kg)/lower background GENIUS detector. Benchmarks G, F, C, A and J can 
be reached by 1 tonne discriminating detectors, such as CryoArray, or ZEPLIN-MAX, discussed above. Some 
of the benchmarks within reach of the direct detection experiements will not be reached by accelerators in this 
decade, illustrating the complementary nature of direct detection results. 


VII. INDIRECT DARK MATTER SEARCHES 


Indirect dark matter signals arise from enhanced pair annihilation rates of dark matter particles trapped in 
the gravitational wells at the centers of astrophysical bodies. The different signals can be classified according 
to the nature of the emitted particles. 

While most of the neutralino annihilation products are quickly absorbed, neutrinos may propagate for long 
distances and be detected near the Earth’s surface through their charged-current conversion to muons. High- 
energy muons produced by neutri nos f r om t he center s of t he Sun and E arth are, therefore, prominent signals 
for indirect dark matter detection |18£, 190, 191, 192, 193, 194, 195, 1961. 
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FIG. 10: Elastic cross sections for spin-independent neutralino-proton scattering. T he p redictions of SS ARD (blue crosses) 
and Ne utdr iver (red circles) are compared. Projected sens itivi ties for CDMS II [18t] or CRESST USm (solid); GE¬ 
NIUS ,187] (dashed); CryoArray \l8(^ or ZEPLIN-MAX [181] (dot-dashed) are also shown. (See Ref. ^.) Further 
theoretical and experimental direct detection data can be plotted using an interactive web plotter ^18^ . A cross section 
pb at a WIMP mass of 100 GeV (in line with current experimental sensitivities) would give an integrated 


of Ax 10 
event rate of 1.0 (T2) 


d for a Ge (Xe) target with a detection threshold of 10 keVr. A cross section of 10 ^ pb 


would give an integrated event rate of 1.0 (1.2) (100 kg) year for Ge (Xe). 
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FIG. 11: Muon fluxes from neutrinos originating from relic annihilations ins ide the Sun. Avo roximate sensitivit ies o f 
near future neutrino telescopes = 1 0^ km“^ yr~^ for AMANDA 11 [191], NESTOR [198], and ANTARES \l9^] , 


and = 1 km yr for IceCube ,20Cj) are also indicated. (From Ref. 


W) 


Muon fluxes for each of the benchmark points are given in Fig. using Neutdriver with a fixed constant 
local density po = 0.3 GeV/cm^ and neutralino velocity dispersion v = 270 km/s (for further det ails, see f|166| ). 
For the points considered, rates from the Sun are far more promising than rates from the Earth [166, 195|| . For 
the Sun, muon fluxes are for the most part anti-correlated with neutralino mass, with two strong exceptions: 
the focus point models (E and F) have anomalously large flrrxes. In these cases, the dark matter’s Higgsino 
content, though still small, is significant, leading to annihilations to gauge boson pairs, hard neutrinos, and 
enhanced detection rates. 

The exact reach of neutrino telescopes depends on the salient features of the particular detector, e.g., physical 
dimensions, muon energy threshold, etc., and the expected characteristics of the signal, e.g., angular dispersion, 
energy spectrum and source (Sun or Earth). Two sensitivities, which are roughly indicative of the potential of 
upcoming neutrino telescope experiments, are given in Fig. For focus point model E, where the neutralino 
is both light and significantly different from pure Bino-like, detection is possible already in the near future at 
AMANDA II 1^, NESTOR and ANTARES [|^. Point F may be within reach of IceCube as 
the neutralino’s significant Higgsino component compensates for its large mass. For point B, and possibly also 
points I, G, C, and L, the neutralino is nearly pure Bino, but is sufficiently light that detection at IceCube may 
also be possible. 
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FIG. 12: The integrated photon flux 3 > 7 (£'th) as a function of photon energy thresho ld E th for photons produced by relic 
annihilations in the galactic center. A moderate halo parameter J = 500 is assumed i204j. Point source flux sensitivities 
for various gamma ray detectors are also shown. (From Ref. ^.) 


As with the centers of the Sun and Earth, the center of the galaxy may attract a significant overabundance 
of relic dark matter particles [ ^01 , 202, 202]. Relic pair annihilation at the galactic center will then produce an 
excess of photons, which may be obs erved in gamma ray detectors. While monoenergetic signals from xx 77 
and XX jZ would be spectacular [204|, they are loop-suppressed and rarely observable. 

Fig. |l^ shows t he i ntegrated photon flux < 1 ) 7 (Eth) in the direction of the galactic center, computed following 
the procedure of |1^|. Estimates for point source flux sensitivities of several gamma ray detectors, both current 
and planned, are also shown. The space-based detectors EGRET, AMS /7 nnd GLAST can detect soft photons, 
but are limited in flux sensitivity by their small effective areas. Ground-based telescopes, such as MAGIC, 
HESS, CANGAROO and VERITAS, are much larger and so sensitive to lower fluxes, but are limited by higher 
energy thresholds. These sensitivities are not strictly valid for observations of the galactic center. Nevertheless, 
they provide rough guidelines for what sensitivities may be expected in comi ng years. For a discussion of these 
estimates, their derivation, and references to the original literature, see [195|. 

Fig. 1^ suggests that space-based detectors offer good prospects for detecting a photon signal, while ground- 
based telescopes have a relatively limited reach. GLAST appears to be particularly promising, with points I and 
L giving observable signals. One should keep in mind that all predicted fluxes scale linearly with J, and for a 
different halo profiles may be enhanced or suppressed by up to two orders of magnitude. Such an enhancement 
may lead to detectable signals in GLAST for almost all points, and at MAGIC for the majority of benchmark 
points. 

Relic neutralino annihilations in the galactic halo may also be detected through positron excesses in space- 
based and balloon experiments |20f;, 206, p0^208|. Ref. |166| also estimated the observability of a positron 
excess, following the procedure advocated in IM For each benchmark spectrum, one first finds the positron 
energy Eopt at which the positron signal to background ratio S/B is maximized, and then requires S/B at 
L'opt to be within reach of the experiment. The sensitivities of a variety of experiments have been estimated 
in 11951. Among these experiments, the most promising is AMS [ ^09| , the anti-matter detector to be placed 
on the International Space Station. AMS will detect unprecedented numbers of positrons in a wide energy 
range. We estimate that a 1% excess in an fairly narrow energy bin, as is characteristic of the neutralino 
signal, will be statistically significant. Unfortunately, all benchmark points yield positron signals below the 
AMS sensitivity |166|. Similar rather pessimistic conclusions were derived in |21C, 211|. Of course, one should 
be aware that as with the photon signal, positron rates are sensitive to the halo model assumed; for clumpy 


halos 1 2 ^ I, the rate may be again enhanced by orders of magnitude 208 


Other possible indirect dark matter signals are antiprotons [193, ^13 , 214] and antideuterium [215], but those 
were not pursued during Snowmass 2001. 


VIII. DISCUSSION AND CONCLUSIONS 

Fig. 1^ shows a compilation of many pre-LHC experiments in astrophysics, as well as particle physics at both 
the energy frontier and lower energies. The signals considered, the projected sensitivities, and the experiments 
likely to achieve them, are discussed in |Q|. On the particle physics side, the signals considered were super- 
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FIG. 13: Estimated reaches of various high-energy collider and low-energy precision searches, direct dark matter searches, 
and indirect dark matter searches before the LHC begins operation, for tan (3 = 10. The projected sensitivities used are 
given in Ref. Q. The darker shaded (green) regions are excluded by the requirement that the LSP be neutral (left) and 
by the LEP chargino mass limit (bottom and right). The regions with potentially interesting values of the LSP relic 
abundance: 0.025 < ^ 1 (light-shaded, yellow) and 0.1 < < 0.3 (medium-shaded, light blue), have also been 

delineated. The regions probed extend the curves toward the forbidden regions. (From Ref. Q.) 


symmetry searches at LEP [216| and the Tevatron 217 , 218, 219, 220, 221|, the improved measurement of the 
B —> Xs'y branching ratio at B-factories, as well as the projected final sensitivity of the Brookhaven — 2 
experiment. On the astrophysics side, the hgure shows the projected reach of the upcoming direct dark matter 
detection experiments, as well as the multitude of other experiments to detect indirect neutrino, photon or 
positron signals from neutralino annihilations, as discussed earlier. 

Several striking features emerge from Fig. |^. First, we see that, within the minimal supergravity framework, 
nearly all of the cosmologically preferred models will be probed by at least one experiment^. In the most natural 
regions, all models in which neutralinos form a significant fraction of dark matter will yield some signal before 
the LHC begins operation. 

Also noteworthy is the complementarity of traditional particle physics searches and indirect dark matter 
searches. Collider searches require, of course, light superpartners. High precision probes at low energy also 
require light superpartners, as the virtual effects of superpartners quickly decouple as they become heavy. 
Thus, the LEP and Tevatron reaches are confined to the lower left-hand corner, as are, to a lesser extent, the 
searches for deviations in i? ^ Xs"f and a^. These bounds, and all others of this type, are easily satished in 
the focus point models with large mg, and indeed this is one of the virtues of these models. However, in the 
focus point models, all of the indirect searches are maximally sensitive, as the dark matter contains a significant 
Higgsino component. Direct dark matter probes share features with both traditional and indirect searches, and 
have sensitivity in both regions. It is only by combining all of these experiments, that the preferred region may 
be completely explored. 

Finally, these results have implications for future colliders. In the cosmologically preferred regions of param¬ 
eter space with 0.1 < < 0.3, all models with charginos or sleptons lighter than 300 GeV will produce 

observable signals in at least one experiment. This is evident for tan/3 = 10 in Fig. |^. In Fig. |^, we vary 
tan/3, fixing M 1/2 to 400 GeV, which roughly corresponds to 300 GeV charginos. We see that the preferred 
region is probed for any choice of tan /3. (For extremely low tan /3 and mg, there appears to be a region that is 
not probed. However, this is excluded by current Higgs mass limits for Ag = 0. These limits might be ev aded 
if Aq is also tuned to some extreme value, but in this case, top squark searches in Run H of the Tevatron |222| 
will provide an additional constraint.) 

These results imply that if any superpartners are to be within reach of a 500 GeV lepton collider, some 
hint of supersymmetry must be seen before the LHG begins collecting data. This conclusion is independent of 
naturalness considerations. While our quantitative analysis is confined to minimal supergravity, we expect this 
result to be valid more generally. For moderate values of tan /3, if the dark matter is made up of neutralinos, they 


^ While this is strictly true for low tan /3, at higher tan (3 some of the nr eferr ed region may escape all probes, but this requires heavy 
superpartners and a significant fin e-tu ning of the electroweak scale [195|. Furthermore, the large tan/3 region is also effectively 
probed via the Bs —> signal |jl54] which was not considered here. 
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FIG. 14: As in Fig. ^ but in the (mo, tan/3) plane for fixed M 1/2 = 400 GeV, Ao — 0, and /r > 0. The regions probed 
are toward the green regions, except for ‘I>®°, where it is between the two contours. The top excluded region is forbidden 
by limits on the GP-odd Higgs scalar mass. (From Ref. [^.) 


must either be light, Bino-like, or a gaugino-Higgsino mixture. If they are light, charginos will be discovered. 
If they are Bino-like, light sfermions are required to mediate their annihilation, and there will be anomalies 
in low energy precision measurements. And if they are a gaugino-Higgsino mixture, at least one indirect dark 
matter search will see a signal. For large tan /3, low energy probes become much more effective and again there 
is sensitivity to all superpartner spectra with light superpartners. Thus it appears, on qualitative grounds, that 
all models in which the scalar masses are not widely separated, and the charginos are not extravagantly heavy, 
will be accessible prior to LHC operation. 

The most sensitive tests for SUSY contributions in B decays before the start of the LHC will come from 
sin 2/3 measurements in charmonium modes and CP-asymmetry measurements of i3 —> Xsj and B K*j by 
BABAR and BELLE. For the latter studies the large SM theoretical uncertainties present in the branching- 
fraction measurements that make an extraction of SUSY contributions difficult are absent. In the LHC era 
improved measurements on sin 2/3 and the CP asymmetry in will be carried out. In addition, measure¬ 
ments of branching-fractions and lepton forward-backward asymmetries of 33 —> modes will provide very 

sensitive tests for uncovering SUSY effects. While hadron-collider experiments will focus on ii~ modes, 

asymmetric P-factory experiments will study B Xst^£~ and B decays. At a super B factory 

the precision can be significantly improved in these modes. While the first studies of a will be performed by 
BABAR and BELLE, precise measurements will come from from BTEV and LHCb. In particular, the hadron- 
collider experiments produce high statistics Bg samples, which allow for precise measurements of 7. Additional 
precise measurements of 7 at the T(4S') can only come from a super B factory. By this time a considerable 
reduction of the theoretical uncertainties in the extraction of |I4b/Ucb|, \Vtd\ and |Vts| from data is also ex¬ 
pected, allowing to perform a sensitive model-independent analysis for extracting a new weak phase and new 
contributions to P°i3° mixing. In order to observe the modes B Xgvv^ B or B ^ K*^vv, which 

bear the least theoretical uncertainties of radiative penguin decays, and study their properties a super B factory 
is a prerequisite. So, for at least the next ten years a vivacious B physics program ensures many high-precision 
measurements of CP asymmetries and radiative penguin decay properties that are suitable for extracting SUSY 
contributions. 
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